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We consider a large condensate in a rotating anisotropic harmonic trap. Using the method of 
matched asymptotic expansions, we derive the velocity of an element of vortex line as a function 
of the local gradient of the trap potential, the line curvature and the angular velocity of the trap 
rotation. This velocity yields small-amplitude normal modes of the vortex for 2D and 3D conden- 
sates. For an axisymmetric trap, the motion of the vortex line is a superposition of plane-polarized 
standing-wave modes. In a 2D condensate, the planar normal modes are degenerate, and their su- 
perposition can result in helical traveling waves, which differs from a 3D condensate. Including the 
effects of trap rotation allows us to find the angular velocity that makes the vortex locally stable. For 
a cigar-shape condensate, the vortex curvature makes a significant contribution to the frequency of 
the lowest unstable normal mode; furthermore, additional modes with negative frequencies appear. 
As a result, it is considerably more difficult to stabilize a central vortex in a cigar-shape condensate 
than in a disc-shape one. Normal modes with imaginary frequencies can occur for a nonaxisym- 
metric condensate (in both 2D and 3D). In connection with recent JILA experiments, we consider 
the motion of a straight vortex line in a slightly nonspherical condensate. The vortex line changes 
its orientation in space at the rate proportional to the degree of trap anisotropy and can exhibit 
periodic recurrences. 



I. INTRODUCTION 

The experimental achievement of Bose-Einstein condensation in confined alkaU-atom gases has stimulated great 
interest in the generation and observation of vortices in such systems . Rotating a totally anisotropic harmonic 
trap at an angular frequency f2 can, in principle, generate vortices; they are energetically stable for > flc Jz] . 
There are several other ideas to create vortices in a trapped Bose-Einstein condensate (BEC) ||5|,pd] pO{ . Vortex 
formation in BEC was recently observed experimentally [pT|-p5[ . 

In general, a vortex line in a trapped Bose-Einstein condensate is nonstationary. The vortex line can move as a 
whole, undergo deformation of its shape or perform oscillatory motion like helical waves p6| , p7| . An extensive literature 
exists on vortex dynamics in superfluids |2^. The nonlinear Schrodinger equation (Gross-Pitaevskii model) p9|-p^] 
has served to study the dynamics and reconnection of vortices, their time evolution, and scattering interactions of 
superfluid vortex rings. Vortex precession in a nonuniform light beam has recently been observed and discussed in 
terms of the nonlinear Schrodinger equation p^ . 

The dynamics of a vortex line in a spatially inhomogeneous two-dimension (2D) condensate was considered in p5| , ^ , 
while the problem of curvature-driven motion of a vortex line in a homogeneous superfluid in three dimensions (3D) 
was studied in |]37|| . A normal mode with negative frequency that corresponds to a vortex precession was found 
numerically and analytically for a large 3D disk-shape BEC ^ , and for a small BEC |^ . The motion of vortex 
lines and rings in Bose-Einstein condensates in harmonic traps was studied in 2D and 3D by numerical solution of the 
Gross-Pitaevskii equation |40[ . Minimum energy configurations of vortices in a rotating trap were considered in . 

In a nonrotating trap, the vortex state has a higher energy than the ground-state Bose condensate, so that the 
vortex is thermodynamically unstable However, the vortex (with unit circulation quantum) is dynamically stable 
and can decay only in the presence of dissipation. Dissipative dynamics and the decay time of the vortex state (due to 
the interaction of the vortex with the thermal cloud) in a trapped Bose-condensed gas are discussed in ^] , where the 
friction coefficient is found to be proportional to the temperature. At temperatures relevant to current experiments, 
one can neglect dissipation in studying the normal modes of the vortex because the vortex decay rate is much smaller 
than the frequencies of the normal modes. 

In this paper we consider the dynamics of a vortex line in a zero-temperature condensate in the Thomas-Fermi 
(TF) limit, when the vortex core radius ^ ^ cP /R is small compared to the mean oscillator length d and the mean 
dimension R of the condensate [here, d — [d^dydzY^^ with di = yJh/Muji and trap frequencies uji (i = x,y,z)]. We 
derive a general nonlinear equation for the motion of the vortex that includes the eifects of the trap potential, the 
vortex curvature and the angular velocity of the trap rotation [see Eq. ( |38| ) below]. Linearization of this equation 
around stationary configurations gives rise to the equation for the normal modes of the vortex line. We investigate 
normal modes of the vortex in 2D and 3D condensates. For a 2D condensate there are solutions in the form of helical 
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waves. For nonrotating trap some of the solutions have negative eigenfrequencies (these modes are formally unstable); 
furthermore, in a nonaxisymmetric trap, some solutions can have imaginary eigenfrequencies, implying that a straight 
central vortex line is unstable with respect to finite self-induced curvature. 

In a 3D condensate the spectrum of normal modes becomes discrete. For a vortex near the z axis, the number 
of normal modes with negative frequency depends on the aspect ratio R^/ Rz- A vortex in a disc-shape condensate 
{Rz < R±) has only one mode with negative frequency. However, if we change the aspect ratio to a cigar-shape 
condensate with R± < Rz, more modes with negative frequency appear. Thus it is more difficult to stabilize a vortex 
in a cigar-shape condensate rather than in a disc-shape one. 

The plan of the paper is the following. In Sec. II we derive a general equation of vortex dynamics using the method 
of matched asymptotic expansions. In Sees. Ill and IV, we discuss the normal modes of a vortex line for 2D and 
3D condensates. In Sec. V we investigate normal modes with imaginary frequencies that appear for a vortex in a 
nonaxisymmetric condensate. In the last section we study the motion of a straight vortex line in a slightly nonspherical 
trap. 



II. GENERAL EQUATION OF THE VORTEX DYNAMICS 

Consider a condensate in a nonaxisymmetric trap that rotates with an angular velocity Jl. At zero temperature 
in a frame rotating with the angular velocity fl, the trap potential Vtr is time-independent, and the evolution of the 
condensate wave function 5" is described by the time-dependent Gross-Pitaevskii (GP) equation: 

(^-|^v2 + Vt, + g\^\^ - fi{n) + ihn ■(rxv)^^ = ih^, (i) 

where T4r = [uj'^x^ + (^yU"^ + is the external trap potential, g = inh^a/M > is the effective interparticle 

interaction strength, and is the chemical potential in the rotating frame. 

We assume that the condensate contains a g-fold quantized vortex with the position vector rQ(z,t). In this section 
we use the method of matched asymptotic expansions to determine the vortex velocity as a function of the local 
gradient of the trap potential Vtr, the vortex curvature k and the angular velocity $7, generalizing the two-dimensional 
results obtained by Rubinstein and Pismen |^,^ to the case of a three-dimensional rotating potential. The method 
applies when the external potential does not change significantly on distances comparable with the core size |(7|^ ^ Rj_ 
(this is the TF limit) and when the curvature is not too large {k <^ l/|g|^); it matches the outer asymptotic form of 
the solution of Eq. (|l|) in the vortex-core region (|p — Pq| < \q\£_) with the short-distance behavior of the solution in 
the region far from the vortex core (|/3 — Pol^klO- 

To find the solution in the vortex-core region, one may consider Eq. (|l|) in a local coordinate frame centered at the 
point Tq of the vortex line that moves with the vortex velocity V. In the general case, the vortex line has a curvature 
k that depends on the specific element in question. We introduce a local coordinate system (x, y, z), so that the x axis 
is directed along the vortex normal, the y axis is along the binormal b and z axis is along the tangent t (see Fig. 1). 

vortex line axis of trap rotation 




Fig. 1. Local coordinate system associated with the vortex line. 
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The solution is assumed to be stationary in the comoving frame and satisfies the equation (in the local coordinates): 



= -ihY ■ , (2) 

where the term —kdx arises from the transformation to local coordinates. 

One can remove Jl from this equation by a shift V V— (Jl x tq). In the vortex core region we may seek a solution 
in the form of an expansion in the small parameters C/R± and k^: 

vf *o(p) + *i = [1*0(^)1 -x(p,^)cos^]e'^^-'''(''^-)-'^^ (3) 

where is the condensate wave function with Vu replaced by Vtr(T"o); it satisfies a zero-order equation 

-|^V2 + Vt,.{ro) + .g|v|/oP - ^(17)) *o = , (4) 

and Xj V characterize the perturbation in the absolute value and phase. Physically 'i'o is the analogous wave function 
for a laterally unbounded condensate with chemical potential — Vtr(po)- The polar angle i/i is measured from the 
direction of the vortex normal (n || £) and p is the radial cylindrical coordinate in the local frame. 
The perturbation 5*1 obeys the following equation 

L *D = • V*o + ^*oP • ViKr(ro) + kdx'ifo, (5) 

II n 

where 

L ^ V^*, + ^ [(^(f7) - V^ro) - 2g|*oP) *i - 9^1^] (6) 

n 

is a self-conjugate operator and is the gradient operator in a plane perpendicular to the vortex line. This equation 
is linear in 'J'l and in V; it contains V^Vtr and k as independent sources, so that the velocity of the vortex line is 
a sum of independent contributions due to V^Vtr and k. Also, the function '^o depends only on the coordinates 
in the direction perpendicular to the vortex line; therefore, in the dot product V • V^'o only the component of the 
velocity perpendicular to the vortex line is relevant. We also assume that V has no component along the line. For 
simplicity, one can assume that Vj_Vtr lies along n and derive the vortex velocity as a sum of two independent 
contributions. The final result written in vector form remains valid for arbitrary directions of V^T^r and n. Under 
this assumption, we have V • V^o — Vdy'i'o = V (sin (j)dp + p^^cos^c?^) in polar coordinates. Then, writing 
^I'l — —(xcoscf) + 177! 5*0 1 sin0)e"''^ in terms of the small perturbations x and rj, Eq. has the form: 

dip + -dp + C) X + ^ {m ~ ^tr(ro) - 3.g|*oP) X - - ^l*o|^ 

= ^l*o| (^^-pIV^FtrMl) -A:9,|vl/„| (7) 



%p + -dp + C - 4 ) + ( 5p|*o|9pry + oia.r; - \-, 

p p J ^0 V p 



2Md^^k_, 



We can remove V from these equations with the following gauge transformation: 



M 

■q = fi- —pV (9) 
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Further, for large distances \q\( ^ p ^ R±, we can use gl^foP ~ sI^tfI"^ ~ ~ ^r('"o) and rewrite Eqs. (0) and 
(0) as follows 



2g\^TF\x = p\V±Vt,{po)\ + ^]^^^M^ol (10) 



equivalently, 



X=,^Zr^\'^Mro)\ + TT7Zn^^dp\^o\, (12) 



r 2^2 ^ o gp|V±t4r(r-o)| , qkh^dpl^ol 

{pO^^ + pdp-l)^-^^^-^^kqp+^^^-^. (13) 

In Eq. (|l^), we can omit the last term, which is smaller with respect to the term kqp by the factor C^/p^. As a 
result, for p 3> \q\S. the perturbations have the following asymptotic form: 



|V±l^tr(ro)| , fcn^ ^ , 

In terms of the phase S", the solution (|lj) (the inner expansion in the coordinate frame centered at the vortex line) 
has the form: 

S-,^-lC^^^ + k)MAp)y+^iV + n.r,).r (16) 

The parameters A and V must be determined by matching the solution ( p^ ) with that far from the vortex core. 

To the lowest order in the small parameter £,/R±, Eq. (|^) far from the vortex core reduces to an equation for the 
condensate phase only 

M 

l^-TFpV^S- + V|*TFp • VS - —n- {r X V) I'J'TFp 0, (17) 
where ^' = |^'|e''^. In the frame rotating with the trap and for fl = Qz, the phase has the form 

where Sq is independent of |Q. Under a shift of coordinates r Tq + r, we have 

M / 2 \ 

S^So^-[in. ro) + j^^^^^ iVVM X 1.)) . . . (19) 

Comparison of (16|) and (|l9|) allows us to find the contribution to the vortex velocity due to the trap rotation 
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V-Vo+ , , (VV^tr(ro)xri), (20) 



where Vq is the velocity for a nonrotating trap. 

It is next necessary to find the asymptotic form of 5*0 far from the vortex core. This function 5*0 satisfies the 
following equation (in the shifted frame): 
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Introduce a function $ such that 



Sox = -q {<S?y + <^dy\n\-^TF\ 



^0, -g(*x + *9xln|*TF|'), 



(21) 

(22) 
(23) 



where Sqx — dxSo and ^oj, — dySo. This representation satisfies Eq. (plj) autoniaticaUy. In addition, the condition 



2 • V X (Vi5) = S-oyx - So^y ^ • (S-oyi - SoxV) = 27rq6^^'> (p) 
gives an equation for $ containing a point source at the vortex location 

Vi_ ■ [Vi$ + $V_L hi I^TFp] = 27r(5(2) (p) . 



Hence 



or 



We can put V^^ (e 



- In I^Tf 



in|*rFlv2^y^j25|^'TFp, SO that Eq. (13) becomes 



Vi $e 



„ln ItTf 



2g|*TF|- 



In |*Tf 



27r(5(2) (pj e 



In 1*tf| 



(24) 
(25) 
(26) 
(27) 

(28) 



To find the sohition, we rewrite Eq. ( psj ) in the local coordinate frame associated with the vortex line, taking into 
account the effect of curvature: 



(^l-kdx) f$ei"l*-H)_ _Zl!:^$ei"l*TH =27r5(2) (p)ei" 



(29) 



(30) 



The solution of this inhomogcneous equation is (note that an additional solution of the homogeneous equation does 
not satisfy the boundary conditions at large p and should be omitted): 



2g\^TF\^ 4 



(31) 



where Kq is a modified Bessel function [for small x, we have Kq ~ — In (e'-^a;/2) where C — 0.577... is the Euler 
constant]. Further, under the logarithm we may put V^ Vtr /45|^tfP ~ Hence $ (in the local coordinate 

frame centered at the vortex line) has the short-distance form 



fc2 



(32) 



To make the asymptotic matching, we can express the solution (|T^) in the vortex-core region in terms of $ and then 
compare with the formula (p2[). Using the definitions (p2|), ( p3| ) of the function one can show that the expression 
( p^ ) (for = 0) corresponds to the following function $ (in the coordinate frame centered at the vortex line): 



2 V .gl^'TFp 



In (Ap) + -^x . 
hq 



(33) 
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To verify this expression, we use In j^'^Fp ~ — | V^Vti(ro)|/5|^'TFp and dylnl^TFl"^ « in the local coordinate 
frame where {dx — k) <&. Substituting into ([22|) and (ESh, we obtain 



^0. = -^, (34) 

P 2 V gl^TFr J ^ 
Therefore, in the coordinate frame centered at the vortex line, the remaining contribution to the phase is 

•So = W - rr- p5 hfc \n{Ap) pain(p+ —Vopsmcj) (36) 



which coincides with ( Jlq ) (for = 0). Matching ( |33| ) and ( p2[ ) (at p ~ |g|^) gives an expression for the constant A 
(with logarithmic accuracy): 



ln(Ae) = lnJ-i^ + ^, (37) 



where R± is the mean transverse dimension of the condensate, and for the velocity Vq. 

Finally, in general vector form (in the frame rotating with the trap), the vortex velocity is: 



"'(^"^ ^ [-W^ ^ 'V ^H'^'^V ^ ^ j A.V^..(ro) ' (''^ 

where 6 is a unit vector in the direction to the vortex binormal, i is a tangent vector to the vortex line and Aj^ is 
the Laplacian operator in the plane perpendicular to n. This formula is valid for arbitrary directions of the local 
gradient of the trap potential, the normal to the vortex line and fl. Near the condensate boundary the denominator 
of the first term in this formula goes to zero. Therefore, t x VVtr must also vanish near the boundary, implying that 
t is parallel to VVtr; as a result, the vortex line obeys the boundary condition that its axis is perpendicular to the 
boundary. 

III. NORMAL MODES IN TWO DIMENSIONS 

To understand the implications of Eq. (|38|), it is valuable to consider first the case of a 2D condensate with Q = 
and cj^ = (hence no confinement in the z direction). Let the vector pQ{z,t) — {x{z,t),y{z,t)) describe the time- 
dependent position of the vortex line during its motion, and k be the vector of the principal curvature; then kb = txh 
and (P p^/ds^ = k, where s is the length measured along the vortex line. For small displacements of the line from the 
z axis, we have s ^ z^ t ~ z and k k.cP p^j dz"^. Then using 

z X VVtr = -Maj^yx + Mtolxy (39) 

and 

kb = txk^ -x-^+y—, (40) 
we obtain the following coupled differential equations for x{z, t) and y(z, t) 

dx qh ( 2y d^y\ , f, H 4f7^ y ,^ , 

dy qh ( 2x d^x\ I , ,^ 1^ \ 4ri/i x 




dt 2M \Rl dz^ J y ' V Rl 8/ M{uj^+ul)Rl 
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These equations have solutions in the form of hehcal waves 

X = Ex sm{ujt + Kz + (fo), y — Ey cos{ut + kz + ipo), (43) 
with the following dispersion relation between uj and k (under the logarithm we take k « |k|): 




The associated amplitudes obey the relation 




(45) 



where 



AMRlRln ^ ^ 

^ = , " — , (46) 

qhiRl+Rl)ln{\q\^^^+\4 

is a dimensionless rotation speed. 

For a nonaxisymmetric trap (for example, R^ > Ry) the oscillation frequency becomes imaginary if \J {2 — fl)/ Rx < 

\k,\ < \J {2 — Cl)/Ry, in which case the initial orientation of the vortex line along the z axis is unstable with respect 

to the formation of finite curvature. For O > = 2, however, the oscillation frequency is real (and positive) for any 
K. Thus the trap rotation stabilizes a vortex line that initially lies along the z axis. 

For a straight vortex line (k = 0), the frequency is always real. This unstable normal mode has the most negative 
frequency [we choose the sign in (^) that corresponds to positive-norm solution] with 

^InfeV J . (47) 



MRxRyy y\q\u iH^l+^l), 

For q > and $1 = 0, the vortex moves (around the z axis) counterclockwise in the positive sense. With increasing 
rotation frequency Vl of the trap, the vortex velocity (as seen in the rotating frame) decreases towards zero and 
vanishes at = fim, where the metastable angular velocity of trap rotation is given by 

a„^^^^4±^lnf^V (48) 



4m Vkie 

This value corresponds to the angular velocity of trap rotation at which a straight vortex line at the trap center 
first becomes a local minimum of energy (for < f2m, the central position is a local maximum) For > fim the 
apparent motion of the vortex becomes clockwise. For k = 0, this straight vortex follows an elliptic trajectory along 
the line Vtr = const, as expected from the dissipationless character of the GP equation. 

For a uniform condensate {Rx, Ry oo), Eq. ( p[ ) coincides with the well-known dispersion law of small oscillations 
of a straight vortex line (Kelvin modes): 

U^^±^k^H\qm. (49) 

Note that one can represent the helical wave solution (^) as a sum of two plane-wave solutions: 

xi — Ex cos{kz) sm{ujt + (fio), yi ~ Sy cos{kz) cos{u!t + ipo), (50) 

2^2 = sin(fcz) sin(cjt + + 7r/2), y2 ^ eySin{kz) cos{ujt + ipQ + tt/2). (51) 

One can easily see that Eqs. (|50| ) and (|5l|) are indeed solutions of Eqs. ( ^ ) and (|4^) with the same dispersion law 
u! = uj{k) as the helical wave (M). In fact, the general motion of the vortex line can be represented as a combination 
of plane-wave solutions; helical waves are just one of the possible combinations and hence do not represent a different 
set of solutions. Solutions (|50|), ( |5l| ) have different parity, but the same eigenfrequency (in 2D there is degeneracy). 
In 3D case, the plane-wave solutions are not degenerate and, therefore, in 3D it is impossible to construct a simple 
analog of helical waves. The general vortex motion in 3D is a combination of plane waves (plane- wave solutions in 3D 
exist, at least, for an axisymmetric trap) with different numbers of nodes along the symmetry axis and hence different 
frequencies. 
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IV. DYNAMICS OF A VORTEX IN THREE DIMENSIONS 



Let us consider small displacements of the vortex from the z axis and Jl — fiz. The vortex curvature is proportional 
to the vortex displacement, so one can put fc « under the logarithm in (^8|). Further, for small displacements 

i X VVu ~ M {x {wlzy' - Loly) + y {loIx - ^Izx')) , (52) 

where a prime denotes derivative with respect to z. Then in dimensionless coordinates x RxX, y RyV, z RzZ, 
Eq. (|3^) becomes: 



2MR,Ry\ (l-z2) ; Vkle/ M{ujl+ul)RxRy' 

qh (2{azx'-x) ,\ f R±\ 4f7/i a; 

ax In h-TT - 777^- — (54) 



2MR,Ry\ (1-^2) ; y\q\^J M{Ljl+ujl)RxRy 



where 



are parameters characterizing the trap anisotropy. One can seek solution of these equations in the form 

X = x{z) sin(a;t + ipo), y — y{z) cos{Lut + ipo) 
and obtain the following ordinary differential equations for x{z), y{z) and ui 

2 {(3zy' -y) r, i, , ^ ..t-N 



2 {azx' — x) 



Coy= ^ ' - ax" + fix, (56) 



where we introduce dimensionless angular velocities 



2MR^RyUj ~ AMRlRln 



qn\n[^)' qniRl + Rl)\n(^^) 

A. Stationary configurations 

Consider a nonrotating trap with 51 = 0. In this section, we seek stationary configurations in which the vortex 
line remains at rest (in this case, the contributions to the velocity from the vortex curvature and the trap potential 
compensate each other). To find the stationary configurations we need to solve Eqs. ( |5^ ) and (|5^) with the condition 
uj = 0. The resulting equations for x and y uncouple; for example, the equation for x{z) has the form: 

(1-z^) x" - 2zx' + -a; = 0. (57) 

^ ' a 

The general solution of Eq. (^^ can be expressed in terms of hypergeometric functions, but it is impossible to satisfy 
the boundary conditions that x{z) should be finite at z = ±1 unless 2/a = n{n + I), where n is an integer (n > 0). 
In this case, the solutions reduce to Legendre polynomials 

x(xPn{z). (58) 
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For example, the first three physical solutions are (we ignore n = which corresponds to a = oo): 

xi = Cz, a = I 



(59) 



X2 =e(l-3z2), 



(60) 



X3 — C { z 
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(61) 



If 2/a 7^ n{n + l), the only possible solution of Eq. ( p7\ j is the trivial one with x = 0. The equation for the y coordinate 
has the same solutions y oc Pm{z), if 2//3 = m{m + 1), and y = 0, if 2//3 ^ m{m + 1). That is, if 2/a ^ n{n + 1) and 
2//3 7^ m(m +1), there are no stationary configurations of the vortex line apart from the straight orientation along 
the z axis. 

One should note that the integer n (or to) enumerates the solutions not only for lj — 0; in particular, n represents 
to number of times that the vortex line (precessing with angular velocity a)„) crosses z axis. For an axisymmetric 
trap (a = (3), we can consider a)„ as a function of a; the function a)„ changes sign at a = a„ = 2/n(n + 1). This 
observation allows us to find the number of normal modes with negative frequency at a fixed value of anisotropy 
parameter a. If a > 1 there is only one mode with negative frequency. If i < a < 1 there are 2 such normal modes. 
If i < a < i there are 3 modes, and so on. If a„ < a < a„_i there are n normal modes with negative frequency. 

B. Dynamics of a vortex in a disk-shape condensate Rz <^ R±: Investigation of unstable mode 



In the limit a, /3 3> 1 the approximate solution of Eqs. (55) and ( |56| ) that corresponds to the unstable mode is 

.v2 ■ 



(62) 



2/ = £ 1 



with the corresponding eigenvalue 



2/9/ ' 



= n-3- — {- + - 



-(- - 



(63) 



(64) 



For = the excitation energy is negative and hence formally unstable. If the trap rotates, the solution (|64| ) becomes 
stable at \VL\ > f^m, where 



(65) 



This expression generalizes that for the angular velocity at which a straight vortex at the center of a thin disk-shape 
condensate becomes metastable H, including the corrections of order and I3~'^ . 



C. Dynamics of a vortex in a cigar-shape condensate Rz 3> i?x: Investigation of unstable modes 



In the opposite limit a, /3 ^ 1, the unstable-mode solution of Eqs. (^ and ( |56| ) corresponds to exponential growth 
of the vortex displacement as a function of z. Such a solution is possible in 3D because the condensate is bounded 
along the z axis. 

For simplicity, consider an axisymmetric trap, so that a — (3. In this case, we have only one equation because one 
can seek a solution in the form x{z) = y{z): 

2 iazx' — x) /, ^ 

ujx^—, ^ - ax" + nx, 66) 

(1 - z2) 
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In the limit a ^ 1 Eq. (|66|) has the following approximate solution 



(67) 



UJ 



a 



(68) 



0, where a thin boundary layer appears. The actual 
in a small vicinity of z = and represents a smooth crossover from the region z < into 

in an elongated cigar-shape condensate 



The approximate solution (67) has a nonanalytic behavior at z 
solution differs from (^ 



the region z > 0. Equation (68) yields the metastable angular velocity fi,, 



8a. i?i Vkk 



(69) 



In contrast to Eq. (|65|) for a flattened condensate, this expression becomes very large for highly elongated trap. 
Consequently, it is significantly more difficult to stabilize a vortex in a cigar-shape condensate than in one with a disk 
shape. 



D. Numerical results for 3D 

We have used Eq. ( |66| ) to evaluate the eigenvalues and eigenfunctions for an axisymmetric trap {a — (3). A finite 
trap rotation produces only a shift of eigenvalues by f), so we can set 57 = 0. In addition, solutions of Eq. ( |66| ) can be 
classified as even or odd functions of z. One can enumerate the solutions by the number m of times that the vortex 
line crosses z axis, m = 0, 1, 2, . . .. The lowest (most negative) eigenvalue corresponds to m = 0. 

In Fig. 2 we plot the angular velocity of the vortex precession a) as a function of the trap anisotropy a — R^/R^ for 
TO = 0, 1,2. In appropriate limits, the numerical solution ljq coincides with those found analytically: a)o ~ ~3 — 
for a > 1, and ujq w —1/a for a <^ 1. The next two solutions uji and 0)2 are proportional to a for large a and diverge 
like —1/a for small a. If a > 1 only one mode has a negative frequency (namely Cjo). If | < a < 1 there are 2 such 
normal modes; if ^ < a < ^ there are 3 modes, etc. (these numerical results coincide with those found analytically). 

The more elongated the trap, the larger the number of modes with negative frequencies (see also [Q). This 
conclusion represents one of our main findings. For a disk-shape condensate, the angular velocity f2„j for the onset 
of metastability is smaller than the thermodynamic critical angular velocity fie, with flm — f The situation is 

completely different for a cigar-shape condensate with a — R^^/R^ < 0.26, because 51,„ then becomes larger than il^. 
For comparison, Fig. 2 includes the line uj — —Clc = —5. 

In Fig. 3 we plot the shape of the vortex line for the lowest (unstable) normal mode (to = 0) for different values of 
trap anisotropy a. The function x{z) is an even function of z without nodes and e — x{z = 0). In Fig. 4 and 5 we 
plot the shape of the vortex line for normal modes with one and two nodes for different values of trap anisotropy a. 
In Fig. 4 Xmax = \x{z ^ Rz)\. 
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Fig. 2. Dimentionless frequencies Cj of the first three normal modes of the vortex as a function of the trap anisotropy 
a — R\/ B?^. The lower horizontal line represents the dimensionless thermodynamic critical angular velocity. 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 
0.1 x/e 

Fig. 3. Shape of the vortex line for the normal mode with the lowest frequency (the most unstable mode) for 

different values of the trap anisotropy a. 



11 




-1.0 -0.5 



0.0 

x/x, 



0.5 1.0 



max 



Fig. 4. Shape of the vortex hne for the first odd normal mode for different a. 




-0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 
0.1 x/e 

Fig. 5. Shape of the vortex Une for the second even normal mode for different a. 
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V. NORMAL MODES OF A VORTEX WITH IMAGINARY FREQUENCIES FOR A 
NONAXISYMMETRIC 3D CONDENSATE 



For an axisymmetric trap all small-amplitude normal modes have real frequencies. For a nonaxisymmetric trap 



(a 7^ P), however, Eqs. (55), and (|5^) can have solutions with imaginary frequencies (in certain regimes of trap 
anisotropy). As an example, we recall that a spherical trap has a special normal mode with m = 1 and Q = 0. Let us 
now consider a nearly spherical trap (this geometry is relevant to the JILA experiments ||4^) 

|a-l|<l, |/3-l|<l. 



One can rewrite Eqs. (pq) and (56) as follows 



where 



'^'-'^ry]-HJ;Uv(;), (70) 



Ho = -{2 + d,[{l~z^)d,]}' ° ^ 



10/' 



f3-l \ ^ 2no/^01 



and is a small perturbation. 

The unperturbed equation corresponds to the equation for the normal modes of a vortex in a spherical nonrotating 
trap and, therefore, all eigenfrequencies of the unperturbed equation are real. The eigenvalue cj = of the unperturbed 
equation is degenerate, and there are two solutions that correspond to (2; = 0: 



(71) 



The eigenfunctions of Hq are Legendre polynomials, and these eigenfunctions form a complete basis. Therefore, one 
can apply the usual perturbation theory to solve Eq. (|7^). The matrix elements are given by: 

Vn - V22 - 0, 



v.2 = lh-i + ln], 



To first order in the perturbation, the eigenfrequencies have the form 



5 / V 5 



If, for example, a < I < (3 (namely Rx < Rz < Ry), then the solutions have imaginary frequencies for Q < 5(1 — a): 



Cb = ±i-f, 7 = 5^/( l-a-ifij f/S-l + ijl) >0. (73) 
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Here, the imaginary frequency means a vortex line oriented along the z axis (which is the intermediate principal axis 
in our example) corresponds to an unstable equilibrium. In contrast, a vortex line oriented along the other principal 
axes is stable. If the angular velocity of the trap rotation increases, however, then the solution ( [72| ) becomes real for 
> 5(1 — a) and the vortex line along the z axis becomes stable because of the rotation. 

It is straightforward to consider a general anisotropic trap (not necessarily close to a spherical shape) . The result 
is the following: if the parameters a and /? satisfy the inequality: 

2 

a < < f3, (74) 

n(n + 1) 

where n is a non-negative integer, then there is a normal-mode solution with imaginary frequency that corresponds 
to the nth Legendre polynomial. Moreover, if 

n[n + 1) mym + 1) 

then there are (n — jn + 1) solutions with imaginary frequencies. Increasing the external trap rotation sequentially 
eliminates such solutions. 



VI. MOTION OF A STRAIGHT VORTEX LINE IN A SLIGHTLY NONSPHERICAL 3D CONDENSATE 

In the previous sections, we studied the motion of the vortex line for small displacements of the vortex from 
equilibrium position (normal modes). In this section we solve the general nonlinear equation of the vortex dynamics 
( ^ ) for a slightly nonsphcrical trap. This problem is directly related to a recent JILA experiment involving the 
evolution of an initially straight vortex line in a nearly spherical condensate ]45| ]. In practice, the trap slightly 
deviates from the spherical shape {Rx ^ Ry ^ Rz)- 

For a strictly spherical trap, Eq. ( psf ) has a solution representing a motionless straight vortex line {i \\ VVtr ) that 
passes through the center of the trap, with the shape 

X = ixS, y = jyS, z = 7^s. (76) 

Here s is the length measured along the vortex line starting from the trap center and 7a;, 7y, 7z are the direction 
cosines of the angles between the vortex line and the principal axes a;, y, z [so that 7^ -f 7^ + 7z = 1j ^ = (7a:; 7i/, 7z)]- 
For a slightly anisotropic trap, however, the solution has approximately the same form as ([76|), but the coefficients 7^:, 
7y! 7z become time-dependent. To find a solution to first order in the trap anisotropy, one can omit the curvature of 
the vortex line and put fc = in (|3^). Also we should use the vortex- free condensate density jvPrFp and take to 
be equal to the value for a spherical trap. Using the standard perturbation theory, we obtain the following equations 

, (77) 

-1)1 



(78) 



(79) 



(80) 

Equations of this type are known in classical mechanics as Euler's equations of rigid-body motion; they describe the 
evolution of the angular velocity of free motion of a body with different principal moments of inertia as seen in the 
body- fixed frame Equations ([78|)-(p0|) have the following integrals of motion: 



for the coefficients 7^:, 73^, 72 (here, we assume that there is no trap rotation): 



/7x- 

or, evaluating the integrals, 



?)qn ( R . ^ 2^ 



^qfi , ( R \ , 2 2n 
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7^+7^ + 7' = 1, 



(81) 



2 2 I 2 2 I 2 2 



const. 



(82) 



That is, the ends of the straight vortex move along trajectories that correspond to the intersection of a sphere and 
an eUipsoid with principal axes proportional to Rx, Ry, Rz- In fact, Eq. ( ^2| ) is the equation of energy conservation 
during the vortex motion. 

In the particular case of an axisymmetric trap (for example, ), Eqs. (^-(|^ have the following 

solution 



Iz = 7z(0) = const, 



(83) 



Ix = 7x{0) cos(tji) + 73^(0) sin(tjt). 



(84) 



where 



Ay. 



ly = 7y(0) cos{Lot) - 7a; (0) sm.(ujt), 



■7.(0) In 



R 

w 



bqh I 1 

2M Vm. 



R 



7.(0) In ^ ^ 



(85) 



(86) 



and [7a;(0), 7j,(0), 7^(0) ] fixes the initial orientation of the vortex line. The line precesses around the z axis (the axis 
of symmetry) at a fixed angle of inclination with respect to the z axis. The frequency of this precession depends on 
the inclination and vanishes if the vortex line is perpendicular to the z axis with 7z(0) = 0. 

We now consider a general nonaxisymmetric trap. To be specific, we assume that uix > ujy > (jJz and introduce new 
scaled functions 



In ( v4r 1 ^ [ujI - ujI){ujI - ujI) 7^ 



4a* 



(87) 



- hqh ( R 
0,1 — m 



4/x \\q\i 



hqh, f R 
m 



4m \m 



Then one can rewrite Eqs. (|7q)-(80) as follows 



\J{^1-^1){^1 -^z) Iz 



^x = ~^y^z, 



(89) 



(90) 



(91) 



(5z = -ftx&y 



(92) 



These equations have the following property: if bx, 6y, 5z is a solution of these equations, then if we change sign of 
any two functions (for example, 6x — > —5x, Sy —by, Sz Sz), we obtain another solution. This property can serve 
to construct solutions that satisfy specific initial conditions. Equations (p0|)-(p2|) have three stationary solutions. Two 
of them (the vortex line parallel to the a; or z axis) correspond to a stable equilibrium, while the third solution (the 
vortex parallel to the y axis) is an unstable equilibrium. 

If |(5a;(0)| < 1(5^(0)1, the vortex line oscillates around z axis (this is one of the equilibrium orientations), and one can 
express the solution of Eqs. (p0)-(p3) in terms of Jacobian elliptic functions as follows: 



Sx = 



^Sm + S^iO) cn (^(52(0) + <52(0)t + C, k) 



(93) 
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5y = ±^6m + <52(0) sn [^5m + <52(o)t + C, k) , (94) 



<5, = ±y^<52(o) + 52 (0) dn (y^<52(0) + ^2(0) t + C,k), (95) 

where the modulus k = y^i5^(0) + (5^(0)/ y^^f(0)"+"5^(0) is less than one, and C is a constant that must be chosen to 
satisfy the initial conditions. The solution is a periodic function of time with the period 



yi(0)Taj(0) J" -/T^k^^h?^ ^«(o) + i2(0) 

where K{k) is the complete elliptical integral of the first kind. The solution (p3[)-(p5[) represents a superposition of 
a nonuniform circular motion in a plane perpendicular to the z axis and oscillations along the z axis: Sz oscillates 
within the following segment: 

V6m-sm<\s.\<^5m+siio) (97) 

If |^x(0)| > |^z(0)| (the modulus k is greater than one), the vortex line oscillates around x axis. In this case one 
can use reciprocal modulus transformation (fc • sn(u, k) = sn(A;u, 1/fc), cn(M, k) — dn(fcu, 1/fc), dn(M, k) = cn(ku, 1/A;)) 
and rewrite the solution (p3|)-(|95|) as follows 

= ±^Sm + SliO) dn (^^2(0) + ^2(0) t + C, 1/k) , (98) 
6y = ±^Sm + SliO) sn (^^Sm + 61(0) t + C, 1/fc) , (99) 

5z = + 6l{0) cn (^<52(0) + S^Q) t + C\ 1/k) . (100) 

The solution is a periodic function of time with the period 

+ «|(0) Jo n/<:'-s»iV y'a2(O) + i2(0) 

If |(52;(0)| = |<^z(0)| (fc = 1), the solution reduces to 



S. = ±Sz = , ' ^, (102) 

cosh (^^2(0) + -52(0) t + C) 



Sy = ±^dm + 61(0) tanh (^^61(0) + 51(0) t + c) ; (103) 

during the motion, the vortex remains on a plane through the y axis, oriented along 6x = ±^z (see Fig. 6); it eventually 
lines up along the y axis (which is a direction of unstable equilibrium for the geometry Ux > tOy > uJz that we are 
considering) . 

Finally, one can rewrite these solutions directly in terms of the parameters 7^;, 7j,, that describe the orientation 
of the vortex line. For example, instead of (|9^)-(|9^), we have 



= W 7^(0) + (J_ J) 7g(0) ■ {Lot + C,k), (104) 



ly = iJlliO) + ^2^ 7^(0) • «n (^i + C, fc) , (105) 
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7. = ±\hm + y%3^7y'(0) • dnM + C,k), (106) 

where 



K - [(^^ - '^.')7l (0) + (o;^ - a;2)72 (o)] ' ^ ^ 



— V^^^^ ^(a;2-a;|)7|(0) + (a;i-a;2)72(o) in ^— 




The period of the motion is given by 



-/ 

Jo 



(108) 



^ -i^(fc). (109) 



■\/l — fc2 sin^ 



For a nonaxisymmetric trap we plot trajectories of the end of the vortex Hne in Fig. 6. The plot corresponds to 
the geometry < Ry < Rz- There are two stable orientations (along the x and z axes) and one unstable (along the 
y axis). The shape of the trajectories strongly depends on the initial orientation of the vortex. 




Fig. 6. Typical trajectories of the end of a straight vortex line (that passes through the condensate center) during its 

motion in a slightly nonspherical trap with R^ < Ry < Rz- 
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For an axisymmetric trap, we have lVx = i^y so that fc = 0. Then, using the properties sn(z,0) = sinz, cn(z,0) = 
cosz, dn(z,0) = 1, we can reproduce formulas (|8^)-(86). 

Because the motion is periodic, one can anticipate one or more revivals of the vortex image in the recent JILA 
experiments. The revival time depends on the orientation of the vortex with respect to the symmetry axes and the 
deviations from sphericity. Recently, such revivals of the vortex image (as well as vortex precession ||2^) was seen 
by the JILA group, and their results agree quantitatively with our theory Q. One should note that our analytical 
results for the normal-mode frequencies are valid with the logarithmic accuracy, namely when ln(i?/|(7|^) ^ 1. It is, 
however, straightforward to go beyond logarithmic accuracy and obtain the numerical correction to the logarithm (see 
Ref. 1^). This correction modifies our formulas for the normal-mode frequencies as follows: instead of ln(i?/|q|^), 
one should use ln(i?/|g|f) -I- 0.675 = ln(1.96i?/|q|^). In the JILA experiments, ln(i?/^) « 3.5, and the inclusion of the 
correction improves the quantitative agreement with the experimental observations. 



VII. CONCLUSIONS 



In this paper, we consider the dynamics of a vortex line in 2D and 3D condensates in the TF limit. We took 
into account the nonuniform nature of the system (namely the trap potential), the vortex curvature and a possible 
trap rotation. We derived a general equation of vortex dynamics and investigated various normal modes of the 
vortex line. For an axisymmetric trap, all eigenvalues are real and the motion of the vortex line can be represented 
as a superposition of planar normal modes with different frequencies. In a 2D condensate, the normal modes are 
degenerate, and, as a result, a superposition of planar waves can produce helical waves. In 3D, there is no such 
degeneracy and there is no simple analog of the helical waves. 

An externally applied trap rotation shifts the normal-mode frequencies and makes the vortex locally stable for 
sufficiently large VI. For a cigar-shape condensate, the vortex curvature has a significant effect on the frequency of 
the most unstable normal mode (that with the most negative frequency), and additional modes with (less) negative 
frequencies appear. As a result, it is more difficult to stabilize central vortex in a cigar-shape condensate than in a 
disc-shape one. 

Normal modes with imaginary frequencies can exist for a nonaxisymmetric condensate (both in 2D and 3D). This 
means that the corresponding equilibrium orientation of the vortex is unstable. As an example of the solution of 
the general nonlinear problem of vortex dynamics, we considered the motion of a straight vortex line in a slightly 
nonspherical condensate. The vortex line changes its orientation in space at a rate proportional to the trap anisotropy. 



ACKNOWLEDGMENTS 



This work was supported in part by the National Science Foundation, Grant No. DMR 99-71518, and by Stanford 
University (A.A.S.). We are grateful to B. Anderson, J. Anglin, E. Cornell and D. Feder for valuable correspondence 
and discussions. This work benefitted from our participation in recent workshops at the Lorentz Center, Leiden, The 
Netherlands and at ECT* European Centre for Theoretical Studies in Nuclear Physics and Related Areas, Trento, 
Italy; we thank H. Stoof and S. Stringari for organizing these workshops and for their hospitality. 



[1] M. H. Anderson et ai, Science 269, 198 (1995). 

[2] C. C. Bradley et ai, Phys. Rev. Lett. 75, 1687 (1995). 

[3] K. B. Davis et ai, Phys. Rev. Lett. 75, 3969 (1995). 

[4] B. Jackson, J. F. McCann, and C. S. Adams, Phys. Rev. Lett. 80, 3903 (1998). 
[5] R. Dum et ai, Phys. Rev. Lett. 80, 2972 (1998). 

[6] K. P. Marzlin, W. Zhang, and E. M. Wright, Phys. Rev. Lett. 79, 4728 (1997). 

[7] G. Baym and C. J. Pethick, Phys. Rev. Lett. 76, 6 (1996). 

[8] F. Dalfovo and S. Stringari, Phys. Rev. A 53, 2477 (1996). 

[9] A. A. Svidzinsky and A. L. Fetter, Phys. Rev. Lett. 84, 5919 (2000). 
[10] D. L. Feder, C. W. Clark, and B. I. Schneider, Phys. Rev. Lett. 82, 4956 (1999). 
[11] B. Jackson, J. F. McCann, and C. S. Adams, Phys. Rev. Lett. 80, 3903 (1998). 



18 



E. L. Bolda and D. F. Walls, Phys. Lett. A 246, 32 (1998). 

T. Winiecki, J. F. McCann and C. S. Adams, Phys. Rev. Lett. 82, 5186 (1999). 

J. R. Anglin and W. H. Zurek, Phys. Rev. Lett. 83, 1707 (1999). 

R. J. Marshall et al, Phys. Rev. A 59, 2085 (1999). 

P. D. Drummond and J. F. Corney, Phys. Rev. A 60, R2661 (1999). 

L. Dobrek et al, Phys. Rev. A 60, R3381 (1999). 

B. Jackson, J. F. McCann, and C. S. Ada ms, Phys. Rev. A 60 , 4882 (1999). 
M. Olshanii, and M. Naraschewski, e-print cond-mat/981131^. 



J. Ruostekoski, B. Kneer and W. P. Schleich, e-print :ond-mat/9908095 
M. R. Matthews et al, Phys. Rev. Lett. 83, 2498 (1999). 
K. W. Madison et al, Phys. Rev. Lett. 84, 806 ( 2000). 



K. W. Madison et al, e-print cond-mat/0004037 



F. Chevy, K. W. Madison an d J. Dalibard, e-prin t |2ond-mat/0005221 
B. P. Anderson et al, e-print sond-mat / 000536^ . 



L. P. Pitaevskii, Zh. Eksp. Teor. Fiz. 40, 646 (1961) [Sov. Phys. JETP 13, 451 (1961)]. 

A. A. Svidzinsky and A. L. Fetter, Phys. Rev. A 58, 3168 (1998). 

R. J. Donnelly, Quantized Vortices in Helium II (Cambridge University Press, Cambridge, 1991). 

F. Lund, Phys. Lett. A 159, 245 (1991). 

S. Rica and E. Tirapegui, Physica D 61, 246 (1992). 

S. Rica and E. Tirapegui, Phys. Rev. Lett. 64, 878 (1990). 

J. Koplik and H. Levine, Phys. Rev. Lett. 71, 1375 (1993). 

J. Koplik and H. Levine, Phys. Rev. Lett. 76, 4745 (1996). 

Yu. Kivshar et al, Optics Comm. 152, 198 (1998). 

B. Y. Rubinstein and L. M. Pismen, Physica D 78, 1 (1994). 
E. Lundh, and P. Ao, Phys. Rev. A 61, 063612 (2000). 

L. M. Pismen and J. Rubinstein, Physica D 47, 353 (1991). 
R. J. Dodd et al, Phys. Rev. A 56, 587 (1997). 

M. Linn and A. L. Fetter, Phys. Rev. A 60, 4910 (1999); 61, 063603 (2000). 

B. Jackson, J. F. McCann, and C. S. Adams, Phys. Rev. A 61, 013603 (2000). 

Y. Castin and R. Dum, Eur. Phys. J. D 7, 399 (1999). 

P. O. Fedichev and G. V. Shlyapnikov, Phys. Rev. A 60, R1779 (1999). 

A. A. Svidzinsky and A. L. Fetter, Physica B 284-288, 21 (2000). 

J. J. Garcia- RipoU and V. M. Perez- Garcia, Phys. Rev. A 60, 4864 (1999). 

E. A. Cornell, private communication. 

See, for example, D. Kleppner and R. Kolenkow, An Introduction to Mechanics (McGraw-Hill, New York, 1973), Sec. 7.7. 



19 



